NASA— TM—1 12537 





V/N 

C 

c 


& 


yju-vv- Trt 

, T V ? - ^ 


SATELLITE ANALYSES OF CIRRUS CLOUD PROPERTIES 
DURING THE FIRE PHASE-II CIRRUS INTENSIVE FIELD OBSERVATIONS 

OVER KANSAS 


Patrick Minnis 1 , David F. Young 2 , Patrick W. Heck 2 , Kuo-Nan Liou 3 , and 

Yoshihide Takano 3 


!NASA Langley Research Center, Hampton, VA, USA 23665-5225 
2 Lockheed Engineering and Sciences Co., Hampton, VA, USA 23666 
3 University of Utah, Salt Lake City, UT, USA 82101 


Presented at 


11th International Conference on Clouds and Precipitation 

Montreal, Canada 
August 17-21, 1992 



SATELLITE ANALYSES OF CIRRUS CLOUD PROPERTIES 
DURING THE FIRE PHASE-II CIRRUS INTENSIVE FIELD OBSERVATIONS OVER KANSAS 


Patrick Minnis 1 , David F. Young 2 , Patrick W. Heck 2 , Kuo*Nan Liou^, and Yoshihide Takano^ 


^NASA Langley Research Center, Hampton, VA, USA 23665-5225 
2 Lockhecd Engineering and Sciences Co., Hampton, VA, USA 23666 
^University of Utah, Salt Lake City, UT, USA 82101 


L INTRODUCTION 

The First ISCCP (International Satellite Cloud 
Climatology Project) Regional Experiment (FIRE) Phase II 
Intensive Field Observations (IFO) were taken over southeastern 
Kansas between November 13 and December 7, 1991, to determine 
cirrus cloud properties. The observations include in situ 
microphysical data; surface, aircraft, and satellite remote 
sensing; and measurements of divergence over meso- and smallcr- 
scale areas using wind profilers. Satellite remote sensing of cloud 
characteristics is an essential aspect for understanding and 
predicting the role of clouds in climate variations. The objectives 
of the satellite cloud analysis during FIRE are to validate cloud 
property retrievals, develop advanced methods for extracting 
cloud information from satellite- measured radiances, and provide 
multiscale cloud data for cloud process studies and for verification 
of cloud generation models. This paper presents the initial results 
of cloud property analyses during FIRE-II using Geostationary 
Operational Environmental Satellite (GOES) data and NOAA 
Advanced Very High Resolution Radiometer (AVHRR) 
radiances. 

2. DATA 

GOES visible (VIS, 0.65 pm) data taken every half hour 
t the nominal 1-km resolution were averaged to obtain an 
^ "effective 4-km resolution to match the corresponding 4-km x 8-km 
infrared window (IR, 11.2 pm) radiances. Each scan line of the IR 
data was duplicated to achieve an effective 4-km resolution. 
Radiances at 12.7 pm (channel 7) were also taken each half hour 
at a 14-km resolution using the GOES multispectral imager (MSI). 
Data from the MSI channel 12 (3.94 pm) were taken with the VIS, 
IR and channel 7 radiances at 20 minutes after the hour at 3-hour 
intervals. The MSI data were duplicated to obtain the same 
resolution as the IR data. AVHRR channels 1 (0.67 pm, VIS), 3 
(3.73 pm), 4 (10.8 pm, IR), and 5 (12.0 pm) were also taken at a 1- 
km resolution from the NOAA-11 and NOAA-12. These satellites 
have ascending nodes at 1430 and 0730 local time (LT), 
respectively. 

In addition to standard 12-hourly soundings at the 
National Weather Service (NWS) stations, temperature and 
humidity profiles were taken from hourly Cross-Chain Loran 
Atmospheric Sounding System (CLASS) soundings at several sites 
in the vicinity of Coffeyville, KS, the hub of the experiment (Fig. 
1). Cloud-top heights were derived from lidar backscattering 
profiles. Surface lidars were located at the hub; Parsons, KS; and 
in northeastern Oklahoma. A down-looking lidar was flown on 
the NASA ER-2 at various times during FIRE-II. 

3. METHODOLOGY 

a. Visible and infrared retrievals 

During the daytime, cloud amount and height are derived 
'ith several techniques. One method uses the VIS and IR 
lannels in a manner similar to that used to analyze the ISCCP 

.data. The VIS reflectance, p v , is interpreted with a 

parameterization of radiative transfer calculations using the 
scattering characteristics of a specified particle size distribution. 


Only two distributions are considered here, a water-droplet 
distribution having an effective radius of 10 pm and an effective 
variance of 0.05, and a cirrostratus (CS) hexagonal ice-crystal 
distribution (Takano and Liou, 1989). The water-droplet (WD) 
distribution, close to that used by the ISCCP, was used to compute 
a Mie-scattering phase function. The CS phase function was 
derived from ray tracing. An optical depth, ty, is derived with 
the parameterization from a given reflectance observation using 
one of the two models. In addition to the cloud reflectance, the 
parameterization accounts for the surface reflectance, Rayleigh 
scattering, and ozone absorption. The VIS optical depth is then 
used to compute IR emissivity of the cirrus cloud which is then 
used to correct the observed equivalent blackbody temperature. 
Tip, to obtain a better estimate of the true radiating temperature 
of the cloud, T c . The value of T c is then converted to cloud height, 
z c , using the nearest available sounding. Estimates of the clear- 
sky temperature, T s , and dear-sky reflectance, pys, are needed to 
perform this analysis. The clear-sky reflectance map was 
developed from minimum reflectances taken during November 
1986. Clear-sky temperature was derived as described by Minnis 
et al. (1987). Details of the technique, the parameterizations, and 
the models are given by Minnis (1991). 

For intercomparison with the surface and aircraft data, 
only data within a small area around the sensor are extracted 
from the image. To provide a larger scale cloud-parameter 
dataset for model intercomparison, the data are analyzed half- 
hourly on a 0.5° grid covering the area between 32°N and 42°N and 
between 93°W and 103°W. The VIS and IR data are converted to 
two-dimensional histograms following the methods of Minnis et 
al. (1990). 

b. Multispectral infrared retrievals 


Another technique applicable both day and night uses 



Fig. 1. Locations of Coffeyville (O), Parsons (P), University of 
Wisconsin lidar (L), NWS sites (O), and CLASS sites (A). 



brightness temperature differences that arise from variations of 
cloud optical properties with wavelength and the nonlinearity of 
the Planck function with wavelength. Estimates of cloud 
emissivity, cloud height, and cloud microphysics can be derived 
^eoretically from an analysis of multispectral infrared data HR 
1 channels 7 and 12 (GOES) or channels 3, 4, and 5 (AVHRR)l by 
ung advantage of the observed temperature differences. These 
differences have been examined both theoretically and with 
observations in attempts to derive various cloud properties from 
the satellite data (e.g., Allen et al., 1990; Stone et al., 1990; Parol 
et al., 1991). As in the earlier methods, the technique used here 
attempts to find the best match between the observations and 
theoretical calculations. 

Radiative transfer calculations using an adding doubling 
model were performed for a range of particle sizes and shapes. 
Cirrus clouds were modeled using ice spheres with effective radii, 
r = 4, 8, 16, 32, and 64 pm, and hexagonal crystals having the CS 
and cirrus uncinus (CU) distributions (Takano and Liou, 1989). 
Surface temperatures, Tg, were varied from 260 K to 320 K. Cloud 
temperatures ranged from 210 K to 270 K. Atmospheric 
attenuation of the IR radiances was estimated in the same manner 
used by Minnis (1991). The channel-7 and -12 atmospheric 
absorption optical depths were 2 and 0.75 times the corresponding 
IR values, respectively. The channel-3 and -5 optical depths were 
set at 1.25 and 0.75 times the channel-4 values, respectively. The 
computations were performed for both a wet and dry atmosphere. 

The results were parameterized as 
3 

q = 2 a n ln(t v / p) n + 24 ln(AT sc ) + 35 p, 0) 

n=0 

where ej is the effective emittance for channel l, AT SC = T s - T c , 
and p = cost), and 8 is the satellite zenith angle. Two sets of 
coefficients, aj, were determined for each channel, one for AT SC > 
5K and another for AT SC < 5K. The effective emittance includes 
both the absorption and scattering effects of the cloud. The 
arameterization fits the radiative transfer results with an rms 
rror of + 0.02 and ± 0.03 for large and small AT sc ’s, respectively. 
The VIS optical depth is = Qexv / Qexh where Q C x ls the 
extinction efficiency and is the extinction optical depth for 
channel i. 



Fig. 2. Comparison of surface lidar and satellite-derived cloud- 
center heights over Parsons, KS, 


A given radiance observation is modeled as 

Bi(D = (1 - C) BiCTsi) + C[ejBj(T c ) + (1 - Ei)Bi(T si )], (2) 

where Bj is the Planck function at the center of channel 1 and C is 
the cloud fraction. 

Similarly, calculations were performed for solar 
reflectance in the VIS and GOES channel 12 for the CS and CU 
distributions and for Mie-scattering water droplet distributions 
having an effective variance of 0.1 and r = 2, 4, 6, 8, 12, 16, 32, and 
64 pm. The results were used to develop a set of bidirectional 
reflectance models, pi(ki,Ti,8 0 ,0 ,y ), where kj is the 
microphysical model, 0 O is the solar zenith angle, and y is the 
relative azimuth angle. Thus, in the daytime, for channels 12 
(GOES) or 3 (AVHRR), the observed radiance is approximated as 

BiCT) = BiCT) + PoEilPiC + (1 - C)psil, (3) 

where p 0 = cos0 o , Ei is the solar constant at the central 
wavelength of channel i, and p s | is the clear-sky reflectance for 
channel i. Values for p s i are estimated at 0.05 and 0.10 over land 
and water, respectively. The non-unit emittance of these surfaces 
is included in the clear-sky measurements of T s at night. The 
separate modeling of the emission and reflectance for these 
channels accounts for the diffuse surface source and the beam solar 
source. 

4. RESULTS 

a. Preliminary VIS-IR results 

An initial comparison of weighted cioud-center heights 
from the NASA Langley lidar (J. M. Alvarez, personal 
communication) with the cloud heights derived using VIS and IR 
GOES data over a small area surrounding Parsons, KS, is shown In 
Fig. 2. These data were taken at various times during FIRE-II. As 
in Minnis (1991), the CS model produces cloud center heights 
which are within 0.2 km of the lidar-derived values. The WD 
model yields underestimates of the cloud heights. Additional 
data covering more hours and sites are being analyzed. 

The prime case study day for FIRE-II is November 26, 
1991. During this day, thin cirrus replaced clear skies over the 
hub around noon and produced a halo around the Sun. A broken 
cirrus layer had developed around 2000 UTC (Universal Time 
Coordinated) centered near 9 km. The surface lidars also 
indicated a lower deck at 6.5 km with multiple levels reported by 
surface observers to the north and west. Multilayered clouds 
prevailed with a low-level cloud deck moving in from the south. 
Both a large-scale network and an Inner network of frequent (3-6 
hr) rawinsonde launches and extensive wind profiler, aircraft, and 
surface measurements provide an exceptional dataset for defining 
and monitoring the dynamics of this developing cyclone over the 
south-central Great Plains. 

Figures 3 and 4 show sequences of cloud height and optical 
depth, respectively, resulting from the gridded analyses for four 
times during November 26. A clear zone around the hub at 1600 
UTC filled in by 1900 UTC and remained overcast through 2100 
UTC. High clouds predominate over Kansas and Oklahoma at 
1600 UTC followed by a mix of cloud layers which give rise to the 
lower cloud heights during the afternoon. Low clouds (zc ^ 2 km), 
evident over north Texas in the morning, spread northward during 
the day. They are probably obscured by the higher clouds over 
southeastern Kansas. Thin cirrus clouds (t v < 4; z c > 7 km) are seen 
over part of eastern Kansas and much of Oklahoma at 1600 UTC 
while denser, high-level clouds are found over northern and 
central Kansas. By 1900 UTC, the thin cirrus clouds are mainly 
over southeastern Kansas with average cloud top heights over 7 
km. At 2030 UTC, thin cirrus clouds with mean heights over 9 km 
pervade over eastern Kansas, Missouri, and central Arkansas. 
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Fig 5. IR multispectral histogram for data over Coffey ville at 
0920 UTC on December 5, 1991. The heavy lines represent 
model results for effective ice particle radius of 16 pm. 
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